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Abstract. In order to better understand the occurrence of phase transitions, 
we adopt an approach based on the study of energy landscapes: The relation 
between stationary points of the potential energy landscape of a classical many- 
particle system and the analyticity properties of its thermodynamic functions is 
studied for finite as well as infinite systems. For finite systems, each stationary 
point is found to cause a nonanalyticity in the microcanonical entropy, and the 
functional form of this nonanalytic term can be derived explicitly. With increasing 
system size, the order of the nonanalytic term grows unboundedly, leading to 
an increasing differentiability of the entropy. Therefore, in the thermodynamic 
limit, only asymptotically flat stationary points may cause a phase transition to 
take place. For several spin models, these results are illustrated by predicting 
the absence or presence of a phase transition from stationary points and their 
local curvatures in microscopic configuration space. These results establish a 
relationship between properties of energy landscapes and the occurrence of phase 
transitions. Such an approach appears particularly promising for the simultaneous 
study of dynamical and thermodynamical properties, as is of interest for example 
for protein folding or the glass transition. 
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1. Introduction 

Speaking about energy landscapes in physics, people typically refer to the study of 
stationary points (i.e., points of vanishing gradient) of an energy function (like the 
Hamiltonian or the potential) of a many-particle system. This concept has proved 
particularly useful for the investigation of dynamical properties like reaction pathways 
in chemical physics or conformational changes in biophysics: having computed 
local minima of the energy landscape as well as transition states (i.e., stationary 
points with exactly one eigendirection of negative curvature connecting two minima) , 
reaction pathways and reaction rates can be predicted from these data (see ;1. for a 
comprehensive textbook on the subject). 

Since energy landscapes have been proved to be a useful concept in dynamics, it 
is tempting to try to understand also statistical physical properties from the study of 
stationary points of a many-particle energy function: Thermodynamic equilibrium 
describes a physical system in the long-time limit, and therefore a link between 
dynamics and equilibrium statistical physics has to exist by definition. Understanding 
what this link looks like in detail has remained a challenge in theoretical and 
mathematical physics for decades, and the mathematical branch of ergodic theory 
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has emerged from the efforts to tackle this problem. In this paper, we will not dwell 
on these issues, and will only take the mere existence of a link between dynamics 
and equilibrium statistical physics as a motivation: A concept that has been proved 
useful in dynamics might have a good chance of being useful in statistical physics as 
well. Therefore, we will investigate in this paper the effect of stationary points of a 
many-particle energy function on statistical physical properties, in particular on the 
existence or absence of phase transitions. 

What is interesting and promising about this approach is that it will allow us to 
study both, statistical and dynamical properties from stationary points of an energy 
function. This should prove useful in particular for the study of physical phenomena 
like protein folding or the glass transitions where both, dynamical and statistical 
features play an important role. 



2. Microcanonical statistical physics 

Consider a classical Hamiltonian system of N degrees of freedom, characterized by 
a Hamiltonian function H{p,q) on phase space, where p = {pi, . . . ,pn) is the vector 
of momenta and q = {qi, . . . , qn) is the vector of position coordinates. The most 
fundamental ensemble in equilibrium statistical physics is the microcanonical one, 
and it takes as a starting point the microcanonical entropy 

sw(e) = -^lnf^jv(£), (1) 

where 

f]jv (e) = J dpdq 6{H{p, q)-Ne) (2) 

is the density of states. The integration in ^ is over phase space, and e denotes 
the total energy per degree of freedom. In later sections we will look for signals of 
phase transitions in this microcanonical entropy. Since we have defined the entropy 
as a function of the energy, we will be able to observe only those transitions which 
are driven by the energy (or the temperature, which is the conjugate variable). To 
observe, say, a magnetic field-driven transition, we would need to consider the entropy 
as a function of the magnetization to which the field is conjugate, and all the results 
presented below had to be translated accordingly. 



3. Stationary points of H and nonanalyticities of the entropy sn 

A stationary point {ps,qs) of the Hamiltonian function H is defined as a point of 
vanishing gradient 1^ 

dHip,,q,) = 0. (3) 

The effect of a stationary point on statistical physical quantities is particularly 
pronounced in the microcanonical ensemble, as we will see in the following example. 

I For the moment the discussion is restricted to classical (i.e., non-quantum mechanical) systems. 
We will comment in the conclusions on how to extend the concepts to quantum systems. 
§ For such a gradient to exist, p and q are assumed to be continuous variables. This excluded discrete 
models like the classical Ising model from our discussion. 
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Example 1. Consider the Hamiltonian Junction 

H{p, q) - ^p^ + -\q^, p,qe R. (4) 

This double-well has two minima at (p, q) = (0, ±1) with energy e = H{0, ±1) = —1/4 
and a saddle point at {p,q) = (0,0) with energy e — iJ(0,0) = [see Fig. [I] (left)]. 
The density of states fli{e) displays nonanalyticities precisely at these values of e [Fig. 
m (right)]. 




Figure 1. Plot of the Hamiltonian function Q as a three-dimensional plot (left) 
and as a contour plot (middle). The corresponding density of states Qi{e) is 
nonanalytic at the energy e = — 1 /4 of the minima and at the energy e = of the 
saddle point of H (right). 



The relation between stationary points of the Hamihonian function H and the 
occurrence of nonanalyticities that we observe for this simple low-dimensional example 
turns out to hold true in arbitrary dimension and for stationary points of arbitrary 
index ilj] To derive a quantitative result specifying this relation, it is convenient to 
make the following assumptions: Let H he a, Hamiltonian function of the standard 
form 

N 9 

i—l ^ 

where nii is the mass associated with the ith degree of freedom. The potential F is a 
function of all position coordinates, and it may include external potentials, two-body 
interactions, or arbitrary n-body interactions. Since the kinetic energy term in ([5]) is 
quadratic in the momenta, it is harmless as regards the occurrence of nonanalyticities 
in statistical physics Hence we can focus our attention in the following on the 
potential V and study the configurational density of states 

^Niv) = J dqS{Viq)-Nv) (6) 

and the configurational microcanonical entropy 

sn{v) = ^InilNiv). (7) 

II The index i of a stationary point gs equals the number of negative eigenvalues of the Hessian Sjy 
of V at Qb- 

^ Such quadratic terms give rise to Gaussian integrals in the partition function or density of states 
which can be solved, resulting in analytic contributions to the canonical free energy or microcanonical 
entropy. 
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We further impose the condition that be a so-called Morse Junction, requiring that, 
for all stationary points of V , the Hessian S^y of V has a nonvanishing determinant, 

detiiyfe) 7^0. (8) 

One may argue that this is an insignificant restriction, since Morse functions form 
an open dense subset of the space of smooth functions 0^ and are therefore generic. 
This means that, if the potential V of the Hamiltonian system we are interested in 
is not a Morse function, we can transform it into one by adding an arbitrarily small 
perturbation. One important consequence of the Morse property is that all stationary 
points of such a function are isolated which allows us to study the effect of a single 
stationary point on the configurational density of states Such an asymptotic 

analysis has been reported in [3J: 

Theorem 1. Let y : G — > R fee a Morse Junction with a single stationary point oj 
index i in an open region G C M^. Without loss oJ generality, we assume V{qs) = 0. 
Then there exists a polynomial P oJ degree less than N/2 such that at v = the 
conjigurational density oJ states (0) can be written in the Jorm 

Q^(v) = Piv)+ ^ ^^--^^^ +o(i;^/^--) (9) 
v/|dct[i3yfe)]| 

Jor any e > 0. Here O is the Heaviside step Junction, o denotes Landau's little-o 
symbol Jor asymptotic negligibility, and 

r (_l)'/2„(Af-2)/20(^) fori even, 

[ (_l)(JV-*)/2(_^)(iV-2)/2Q(_^) . 

For a proof of this theorem, the density of states is calculated separately below 
and above v = Q. By complex continuation it is possible to subtract both contributions 
and to evaluate the leading order of the difference. A detailed proof of an even stronger 
result (including higher order terms) is given in [3]. 

The content of Theorem [T] can be summarized as follows: 

(i) Every stationary point qs of V gives rise to a nonanalyticity of the configurational 
entropy sn{v) at the corresponding stationary value v = ~ V{qs)/N. 

(ii) The order of this nonanalyticity is [(TV — 3)/2j, i.e., sn{v) is differentiable 
precisely [{N — 3) /2J -times at v — Vg- 

It is interesting to note that the functional form of the nonanalytic term (fTO|) does not 
depend on the precise value of the index i of the stationary point, and the effect of a 
stationary point on the entropy is more or less independent of the index. Remarkably, 
in spite of this independence, several approaches have been reported in the literature 
in which statistical properties of many-particle systems are estimated exclusively from 
the minima of the energy landscape, disregarding all stationary points of higher index 
[HlGj. It is unclear to the author how these two observations fit together. 

Numerical studies as well as heuristic arguments indicate that the number of 
stationary points of a generic potential V increases exponentially with the number of 
degrees of freedom iV [7]. As a consequence, for large (but finite) N we can expect to 
find a very large number of nonanalyticities of sat. Such behaviour is at pronounced 
variance to the properties of canonical thermodynamic functions of finite systems 
which are known to be always analytic. 
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4. Phase transitions in finite systems? 

In agreement with the intuition produced in Example [U we found a relation 
between stationary points of the many-particle potential V and nonanalyticities 
of the configurational entropy s^- But what is the physical significance of these 
nonanalyticities? 

The microcanonical entropy is a so-called thermodynamic function, and the role 
that this function plays for the microcanonical ensemble is equivalent to that of the 
canonical free energy /^v for the canonical ensemble. The canonical free energy is 
known to be an analytic function for all finite numbers N of degrees of freedom 
[5], and it can develop a nonanalyticity only in the thermodynamic limit N oo. 
These nonanalyticities of /at are of utmost importance in statistical physics and 
thermodynamics, and their presence is often taken as the defining property of a phase 
transition. But does it make sense to consider nonanalyticities of the microcanonical 
entropy as phase transitions of finite systems in the microcanonical ensemble? 

I don't think so, the reason being that these nonanalyticities of sn do by no 
means show the remarkable properties which account for the interest of physicists in 
nonanalyticities of the canonical free energy. For typical many-particle systems, the 
nonanalyticities of /at (if they exist at all) are of very low order, meaning that already 
one of the first few derivatives is discontinuous. Such a nonanalyticity corresponds to 
a dramatic change of the physical properties due to collective effects. In contrast, the 
nonanalyticities of the finite-system configurational entropy typically occur at a huge 
number of different energy values, and their order decreases as the number of degrees 
of freedom becomes larger. The effect on the system's physical properties in that case 
is minor, and this is simply not the dramatic physical phenomenon which we would 
like to call a phase transitionf3 

5. Stationary points of V and phase transitions in the large-A^ limit 

Having observed that stationary points of V can be numerous and that the 
corresponding nonanalyticities of the microcanonical configurational entropy become 
weaker with increasing N, it may appear doubtful whether they are related to the 
occurrence of phase transitions in the thermodynamic limit at all. However, for several 
models, analytic calculations of all stationary points of V have been performed, and 
the results give strong indications that a relation between these stationary points and 
phase transitions in the thermodynamic limit does exist. 

Example 2. One of those models for which all stationary points ofV can be computed 
analytically is the mean-field k -trigonometric model, characterized by the potential 



where A > is a coupling constant and the coordinates Qi G [0, 27r) are angular 
variables. The potential describes a k-body interaction where fc £ N, and the model is 
known to undergo a phase transition for k ^ 2 at a transition potential energy v — A, 

+ Note that here it was not our intention to define the finite-system analogue of a phase transition, but 
attempts of that kind have been pubhshed in the literature. The most widely used definition relates a 
first-order transition in a finite system to a convex intruder in the (otherwise concave) microcanonical 
entropy function, which is equivalent to an S-bend of the microcanonical caloric curve [9l 110). 




(11) 
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whereas no phase transition takes place for k = \. In J 11]/ all stationary points of V 
as well as the corresponding indices have been calculated. It is convenient to arrange 
these data in the form 

N 



a(v) = lim — In 



(12) 



where fij (v) is the number of stationary points Qs of V with index j and with a potential 
energy V{qs)/N not exceeding v. A plot of a is shown in Fig.\^ (left). This quantity, 




Figure 2. Plot of the graph of a as defined in Eq. I|12| l (left) and of the 'flatness 
indicator' ji (right) as functions of v/A for the mean-field fc-trigonometric model 
in the thermodynamic limit. See Example [3] for further details on j£. 



containing exclusively information about the stationary points of V , already signals 
the absence or presence of a phase transition, being smooth for k = 1 in the absence 
of a phase transition, and nonanalytic precisely at the phase transition energy w = A 
when k ^ 2. 

Other models for which the behaviour of a already indicates the presence of a 
phase transition are listed in Table 1 of [H]. Further evidence of a relation between 
stationary points of V and the occurrence of phase transitions in the thermodynamic 
limit is provided by a theorem due to Franzosi and Pettini [T3] . Although this theorem 
is originally stated in the language of topology changes of configuration space subsets, 
it can be rephrased in terms of stationary points of v. Here we give only a sloppy 
reformulation of this result: 

Sloppy theorem 1. Let V be the potential of a system with N degrees of freedom 
and short-range interactions. If some interval [a, b] of potential energies per degree 
of freedom remains, for any large enough N, free of stationary values of V , then the 
configurational entropy s(v) does not show a phase transition in this interval. 

Note that a precise formulation of this theorem requires further technical 
conditions on the potential V (see jl3] for details). 

Although stationary points of V are necessary for a phase transition to occur, 
their presence is by no means sufficient. This is evident from Example O where 
stationary points were found to lie densely in an interval of the potential energy 
(per degree of freedom) axis, but a phase transition takes place at the single value 
of w = A only. To better understand the microscopic origin of a phase transition, 
we might now like to ask: Which of the stationary points of V may give rise to a 
nonanalyticity of the configurational entropy s{v) in the thermodynamic limit? An 
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answer to this question was given in Ref. [14j . The idea behind this result is to sum 
up the nonanalytic contributions to the configurational entropy as given in Eqs. © 
and (jlOp for aU stationary points and to perforin the thermodynamic hmit of this sum. 
Then a bound on the magnitude of this sum is derived which can be interpreted in 
the foUowing way: 

Sloppy theorem 2. The sum of the nonanalytic contributions of the stationary points 
of V to the configurational entropy cannot induce a phase transition at the potential 
energy per particle v — Vt if, in a neighborhood of Vt , 

(i) the number of critical points is bounded by exp(CN) for some C > and 
(ii) the stationary points do not become 'asymptotically flat' in the thermodynamic 



'Asymptotically flat' here means that the determinant of the Hessian of V at the 
stationary points goes to zero in the thermodynamic limit in some suitable sense (see 
[14] for details). This result qualifies a subset of all stationary points of V as 'harmless' 
as what regards phase transitions and leaves only the (we hope few!) asymptotically 
flat ones as candidates for being at the origin of a phase transition. To illustrate 
the power of this theorem, we reconsider the mean-field fc-trigonometric model from 
Example [2l 

Example 3. Without going into the details, we will consider a function je{v) which 
can be computed from the stationary points ofV and has the property of being divergent 
whenever stationary points with stationary value v become asymptotically fiat in the 
thermodynamic limit, and finite otherwise (see '41 for details). For the mean-field 
k-trigonometric model ill]) , a plot of this function is shown in Fig. [H (right). In 
agreement with Sloppy theorem ji is finite for k = 1 where no phase transition 
takes place, but shows a divergence at the transition potential energy v — A for k ^ 2. 
Despite the superficial similarity of the graphs of ji and a in Fig. [2 it is worth pointing 
out that, by virtue of Sloppy theorem ji has a predictive power as what regards the 
occurrence of phase transitions which a is lacking. 

With the results of the present section we have established conditions on the 
microscopic level, i.e., local properties in configuration space, which are relevant for 
the occurrence of a phase transition on the macroscopic level. Such conditions we may 
regard as a way of better understanding the origin of a phase transition. 

Note that throughout Sees. [SHS] we have made reference not to the density of 
states or the entropy, but to the configurational density of states or the configurational 
entropy as defined in ^ and ([7]) . The density of states is related to its configurational 
counterpart by means of a convolution (see Eq. (4.3) of Ruelle's textbook [TS]), and this 
relation allows us to translate our results on the configurational entropy into results 
on the entropy. For finite systems, we found in Theorem [T] that the configurational 
entropy will typically display nonanalyticities of order \_{N — 3)/2j . As a consequence 
of the above mentioned convolution, nonanalyticities of the finite-system entropy will 
therefore be of order — 1, still increasing linearly with the system size A^ [16J . 




For a Hamiltonian function H of standard form (O, every stationary point Qs of 
V is in one-to-one correspondence to a stationary point (0, Qs) of H. As a consequence, 
in the thermodynamic limit our sloppy theorem [2] can be taken over verbatim to the 
entropy: The sum of the nonanalytic contributions of the stationary points of H to 
the configurational entropy cannot induce a phase transition at the energy per particle 
e = Et if, in a neighborhood of £t, the number of critical points is bounded by exp(CA^) 



limit. 
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for some C > and the stationary points do not become 'asymptotically flat' in the 
thermodynamic limit. 

6. Conclusions and outlook 

6.1. Conclusions 

In Sees, ini and [51 we have studied nonanalyticities of the microcanonical entropy of 
finite and infinite systems, and in particular their relation to stationary points of the 
potential V . In summary we have seen that, for systems with standard Hamiltonian 
functions (O, the following statements hold true. 

(i) Stationary points of V cause nonanalyticities of order \_{N — 3)/2j in the 
configurational entropy of finite systems, at pronounced variance to the analytic 
behaviour of the canonical free energy density of finite systems. 

(ii) The number of stationary points is believed to increase exponentially with N for 
generic potentials V . 

(iii) Due to the large number of finite-system nonanalyticities, their relation to phase 
transitions is clearly not one-to-one, but has to be somewhat more subtle. 

(iv) In short-range systems, stationary points of V are necessary for a nonanalyticity 
to occur in the entropy, the configurational entropy, or the free energy in the 
thermodynamic limit. 

(v) Stationary points of asymptotically vanishing curvature may cause a phase 
transition in the thermodynamic limit. 

6.2. Quantum outlook 

The concepts presented, based on stationary points of the potential energy function V ^ 
are of purely classical nature. However, a geometric formulation of quantum mechanics 
(|17j: see [18] for an introduction) provides a suitable framework to take over all the 
key features of these concepts to quantum mechanics. In this geometric framework, 
starting from a Hamiltonian operator H on Hilbert space 7i, an energy expectation 
value function h is defined on the quantum phase space T'{Ti.), i.e., on the complex 
projective space corresponding to 7i. The function h can be shown to have a number 
of stationary points which, remarkably, correspond to the eigenstates of the operator 
H. Further elaboration and application of these concepts in a quantum context is 
currently in progress. 

6.3. Outlook on applications 

Up to now we have mostly emphasized the conceptual aspects of an analysis of 
stationary points of the potential energy V: We have gained insights into the 
analyticity properties of the microcanonical entropy and into the origin of phase 
transitions. Additionally, an analysis of the stationary points of V is of interest in 
applications as well: For dynamical properties, as mentioned in the introduction, 
energy landscape methods have already been used extensively. Together with the link 
between stationary points and statistical physical properties presented in this article, 
it might be promising to investigate both, dynamical and statistical properties of a 
given system simultaneously on the basis of stationary points of the potential energy 
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landscape. Such an approach should prove particularly useful for the study of protein 
folding or the glass transition where dynamical as well as statistical features play a 
role. 

Both phenomena mentioned, the glass transition as well as protein folding, 
have already been investigated from an energy landscape perspective: For the glass 
transition, these studies were pioneered by Stillinger and Weber [5] who studied 
minima of the energy landscape, and generalizations to stationary points of higher 
index were reported by Grigera et al. [19 and others. Similarly, minima and other 
stationary points have been used for analyzing protein folding (see [T] for an overview 
and a list of references). For both types of phenomena, the results suggest that 
studying stationary points of the energy landscape can be an efficient tool, but several 
difficulties have been encountered as well. The hope is that the results presented in 
this article will prove useful for identifying the properties or types of stationary points 
which are most relevant whenever entropic effects play a role. 
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